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ANALYSIS OF HEAT TRANSFER FOR A NORMALLY 
IMPINGING LIQUID-METAL SLOT JET 
by Robert Siegel 
Lewis Research Center 

SUMMARY 

Liquid metals are very effective heat-transfer fluids. For critical local cooling 
applications a useful technique may be to use an impinging liquid-metal jet. In this 
report a two-dimensional analysis is performed for a liquid-metal slot jet impinging 
normally against a flat plate that is heated uniformly. Wall temperatures and heat- 
transfer coefficients are obtained as functions of position along the plate for various 
values of the Peclet number. Liquid metals are low Prandtl number fluids, and as a 
consequence a thermal boundary layer will develop much more rapidly than a viscous 
layer. The fluid can then be approximated as inviscid because the viscous layer is of 
minor importance in the thermal- boundary- layer growth. It is also assumed that turbu- 
lence is small so that molecular conduction is dominating the diffusion of heat. The 
solution is obtained by mapping the jet flow region into a potential plane in which the re- 
gion occupies a strip of uniform width. The energy equation is transformed into poten- 
tial coordinates, and an exact solution is obtained in the strip region. Conformal 
mapping is then used to transform the solution into the physical plane to yield the tem- 
perature distribution along the impingement plate . 


INTRODUCTION 

A possible technique for localized cooling is to use a jet impinging against the heat- 
transfer surface. Liquid metals are very effective heat-transfer fluids and hence may 
prove useful in critical cooling applications. 

For a liquid metal the molecular diffusion of heat is much larger than the molecular 
diffusion of momentum (low Prandtl number fluid), and hence for developing flow and 
heat transfer the viscous boundary layers are much thinner than the thermal boundary 
layers. As a result, the assumption is often made that the viscous- layer development 



can be neglected in computing the heat transfer. The velocity field within the jet region 
can then be determined by using the techniques of inviscid free jet theory. For a two- 
dimensional situation in Cartesian coordinates the free streamlines can be found by con- 
formal mapping (refs. 1 and 2). For an axisymmetric jet, numerical techniques are 
used (refs. 3 and 4). The present study is limited to two dimensions in Cartesian coor- 
dinates (a slot jet) so that an analytical solution can be found. 

In the conformal mapping solution the two-dimensional jet flow region is mapped 
into a potential plane where it occupies a simple region, namely, a strip of constant 
height and infinite length. The convective energy equation and boundary conditions can 
be transformed into potential coordinates and in this coordinate system can be solved by 
available analytical techniques. This coordinate transformation of the energy equation 
has been used to determine the heat transfer for bodies in crossflow (ref. 5). The pres- 
ent analysis is a combination of the inviscid free jet analysis and the energy equation 
transformation technique. In the present problem these techniques reduce the analysis 
to obtaining the temperature distribution in a slab of uniform thickness moving across 
a plane of distributed heat sources . The solution can be found by generalizing a case 
treated in reference 6. 

Results are obtained for the wall temperature and heat-transfer coefficient along the 
impingement plate for various values of the Peclet number. Simple approximate rela- 
tions are obtained at the stagnation point. The stagnation results are compared with 
those in reference 7, which treats gases and ordinary liquids by use of a boundary- 
layer type of analysis. 


SYMBOLS 


b half- width of undisturbed jet 

Cp specific heat of fluid 

h local heat-transfer coefficient along plate 

K modified Bessel function of second kind of order n 
n 

k thermal conductivity of fluid 

m integer 

Nu Nusselt number, h2b/k 

Pe Peclet number, | v^ 1 2b /a 

Pr Prandtl number, C u/k 

r 
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q w heat flux specified at wall 
Re jet Reynolds number, | j 2b p/p 

T dimensionless temperature, tk/bq w 

t temperature 

U dimensionless velocity in x direction, u/lv^J 

U dimensionless fluid velocity vector, u/jv^ | 

u fluid velocity vector 

u, v velocities in x and y directions 

X, Y dimensionless coordinates, x/b and y/b 

x, y rectangular coordinates along and normal to plate 

a thermal diffusivity, k/pC 

r 

r gamma function 

tj dummy variable of integration 

p fluid viscosity 

v fluid kinematic viscosity, p/p 

4 coordinate normal to <£ - plane 

p fluid density 

$ dimensionless potential, <p / 1 v^ | b 
(p potential function 

'k dimensionless stream function, ^// 1 v^ | b 
\fy stream function 

Subscripts: 
w at wall 

°° undisturbed fluid condition 


ANALYSIS 

Geometry 

The two-dimensional flow configuration is shown in figure 1. 
disturbed width 2b impinges normally against a flat plate at y = 


A slot jet with an un- 
0. The plate has a 


3 



uniform heat flux q supplied along it. After turning, the flow moves out in the posi- 
tive and negative x directions in a symmetric fashion and for the inviscid-irrotational 
situation considered here reaches an asymptotic thickness b. From symmetry, only 
the first quadrant of the flow need be considered. In figure 1 some of the boundary 
points have been numbered 1 to 5 for convenience in identification. 


Governing Equations 

The fluid considered in this report is a liquid metal, and since liquid metals have 
very low Prandtl numbers, the viscous diffusion is small compared with the molecular 
diffusion of heat. Hence, in a flow where both the viscous and thermal boundary layers 
are developing simultaneously, there will be little diffusion of vorticity into the thermal- 
boundary-layer region. The viscous layer will be relatively thin, and the thermal layer 
will be essentially all in the inviscid region. A flow that is initially irrotational will 
develop little vorticity and hence can be assumed to remain irrotational. It is also 
assumed that the jet Reynolds number is low enough so that turbulent heat diffusion is 
not important in comparison with the molecular diffusion of heat, which is large for a 
liquid metal. These are standard assumptions for liquid-metal heat-transfer analyses 
in the low Reynolds number range (RePr < 50 in a tube or <500 for crossflow over tube 
bundles), as discussed in reference 8, page 300. 

With these assumptions the fluid velocity can be obtained as the gradient of a poten- 
tial 


u = V cp (1) 

where cp is governed by Laplace's equation, 

V 2 ^ = 0 (2) 

The Cauchy- Riemann equations apply so that the velocity components are related to the 
stream function by 


VL= ^L = ML (3a) 

dx dy 

v = i£ = - M (3b) 

3y dx 
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where the stream function also satisfies Laplace’s equation, 

V 2 xp = 0 

The energy equation for the flow is given by 

pCu • Vt = kV 2 t 

r 


( 4 ) 


( 5 ) 


or, when equation (1) is used, by 


pC p V^ • Vt = kV 2 t (6) 

Boundary Conditions 

In connection with equations (2) and (4) governing the flow, the cp and \p must be 
specified along the boimdaries. Along the axis of symmetry 34, the u component is 
zero, and along the x axis 45 the v component is zero. The entire boundary 345 will 
arbitrarily be designated as the zero streamline. These conditions give 


if = o 

0X 


\ x = 0, y ^ 0 


* =qj 


■'s 


= o 

9y 


> 0 < x ^ , y - 0 


* = 0 ) 


(7a) 


(7b) 


Along the free streamline 12 the pressure is constant, since the jet is exposed to an 
external region at constant pressure. Consequently, the velocity magnitude along 12 is 
equal to | v^ | . The magnitude of the stream function along 12 is found by noting that the 
difference between this stream function and the value i/a = 0 on the zero streamline 
must equal the volume flow. Hence, along the free streamline the conditions for (p 
and rp are 



(8a) 
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^ = v oo b x >y on 12 


(8b) 


Across the width of the undisturbed incoming jet the velocity has a uniform value 
-v^. Hence, from equations (3a) and (3b) f 




0X 

i£ = V( 

0X 


> x, y on 23 


( 9 ) 


At large x the velocity is uniform (u = v^ , v = 0), so that 


£& = v 
3y 


■\ 

oo 

► 


3y 


( 10 ) 


The thermal boundary conditions are the following. There is symmetry of the tem- 
perature distribution about the axis 34 so that 


— = 0 x = 0, y ^ 0 (11) 

0X 

It is assumed that along the free streamline there is negligible heat loss so that 

n g • Vt = 0 x, y on 12 (12) 

Along the solid boundary there is an imposed uniform heat flux q w so that 

— = 0 — x — 00 , y = 0 (13) 

0y k 

Along the cross section of the incoming undisturbed jet the fluid is at uniform tempera- 
ture t^: 


t=t co x, y on 23 


(14) 
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Equations and Boundary Conditions in Dimensionless Form 


The preceding equations and boundary conditions can be placed in dimensionless 
form by defining the following variables: 


X =- 


V = bV 


Y = y 

b 




(Voo |b 


v b 

oo 


tk „ l v =J 2b 

T = Pe = 


The flow and energy equations (1), (2), and (6) become 


U = 


The boundary conditions are 


/vO 

= 0 


^ V$ • VT = V^T 


'k = 0 > x = 0, y s. 0 (x, y on 34) 
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33? 

3Y 


= 0 


^ = 0 > 0 £ x s oo } y = o (x, y on 45) 


il = _ I, 

3Y 


( s $\ 2 + / a $\ 2 _ ^ 

V 3X/ \3Y/ 

* = 1 /> x, y on 12 


h_ • VT = 0 
s 


3<f> 

= 0 

3X 


3* 

= 1 

3X 


T = 

T 

OO 


J 

3* 

"N 

= 1 

3Y 


33? 

= 0 

3Y 

J 


> x, y on 23 


x, y on 15 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


Solution of Flow Problem 

The solution for an inviscid two-dimensional jet slot jet striking a plane is available 
in many textbooks as an example of conformal mapping applied to free jet theory. In the 
course of the solution the jet is transformed into a potential plane as shown in figure 2. 

In this plane the flow has a uniform velocity and is from left to right in a two-dimensional 
channel of unit height. The mapping functions relating the potential and physical planes 
can be obtained as a special case of the results in reference 9. 
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The heated plate 45 in figure 1 is along the positive $ axis (4^ = 0) in figure 2. 
The locations along the positive $ axis are of interest here, as it is the heat -transfer 
behavior along the plate that is desired. From the results in reference 9 the corre- 
spondence between $ and X is given by 


* I 2 1 + U 

n 1 - u 2 


x|* =0 =;( 2tan ' lu + ln 7 ± ^) 


> 0 < U < 1 


(24) 


where U is the dimensionless velocity in the physical plane along the plate in the x 
direction. The results in equation (24) are all that will be needed from the flow portion 
of the problem. 


Formulation for Solution of Energy Equation by Use of Potential Plane 


The strip region in figure 2 is a simple geometry, and hence it provides a conve- 
nient region in which to solve the energy equation. The solution can be evaluated along 
the boundary 45 in figure 2, where the quantities of interest are functions of 4>. Equa- 
tion (24) can then be used to find the corresponding X values so that the heat-transfer 
behavior is obtained along the plate in the physical plane. 

The energy equation, equation (18), has the same form as that of equation (16) of 
reference 10, and this reference gives the details of the transformation from X, Y to 
$, 4' coordinates. Using equation (26) of reference 10 yields the energy equation 


8 2 T | 5 2 T 
3$ 2 3* 2 


Pe 0T = o 
2 0<f> 


(25) 


This is the same energy equation as for convection to a- flow with uniform velocity in a 
parallel- plate channel. The flow would be in the 4> direction with the channel width ex- 
tending across the 4' direction. 

The temperature boundary conditions in equations (19) to (22) must also be express- 
ed in terms of $ and ^ . Equation (19) expresses symmetry about the jet axis. This 
symmetry is preserved in the transformed plane so that the condition becomes 
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— = 0 on 34 

3* 


(26) 


To transform equation (20) use the relation that at a fixed X 


3T _ 3T_ 3£ 3T 3$ 
3Y 3* 3Y 3<f> 3Y 


Applying the conditions from equation (20) along 4' = 0 gives 


-1 - 


3T 

3* 


*=0 


3# 

3Y 


*=0 


From equation (3a) 


U = 


3* 

3Y 


Hence, the boundary condition becomes 


3T 1 


3* U($) 


$, 'J' on 45 


(27) 


where from equation (24) 


; 7T$/2 _ 1 


, 1/2 




1/2 


Since the free streamline 12 is a line of constant , the condition equation (21) becomes 


-52. = 0 ¥ on 12 

3* 


(28) 


From equation (22) 


T = T <f>, ^ on 23 


(29) 


These boundary conditions are summarized in figure 2. 
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In the potential plane (fig. 2) the situation is a channel flow with uniform velocity 
(i.e. , a moving slab) and a uniform entering fluid temperature T . Equation (25) is 
the convection equation with the axial conduction term 3 T/3$ included. The channel 
boundary conditions correspond to insulated walls except for a nonuniform heat addition 
along the positive <£ axis . The solution can be obtained by utilizing some results from 
reference 6. The desired result is the temperature distribution along the plate that the 
jet is impinging against. This can also be expressed in terms of a local Nusselt number 
along the plate. The local heat-transfer coefficient is h = q w /(t w - t^ ). Then the local 
Nusselt number based on the jet width 2b is 


Nu(X) = 

k 


q w 2b _ 2 

(t w -t ") k " TJX) - T, 


(30) 


Solution of Energy Equation 

Equation (8) on page 268 of reference 6 gives the temperature on the surface of a 
slab that is moving past a line source of heat. If superposition is used to obtain the ef- 
fect of a distribution of line sources, the solution can be found for the present problem. 
A difficulty is that the result given in reference 6 is valid only for positive distances 
along the surface away from the line source. A relation valid for both positive and neg- 
ative distances is needed to superpose the line sources and obtain a continuous source 
distribution as given by - 1/U(4>) in figure 2. The details of the solution are given in 
appendix A, and the final result for the wall temperature along the positive $ axis 
(segment 45) is given by 
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m=l 


dr) 


tanhI2 


, 1/2 



2 


[■•raj 

1/2 



s -(Pe/2 )(n-*) 

: di) 

Kf 2 


$ r 

(Pe/4)($-?j) il-|l+(4m7r/Pe) 2 J 


l/2l 


tanhM 


, 1/2 


drj 


■(P e/4)(rf-<f>) ■il+|l+(4m^/Pe) 2 J 


1 / 2 ) 


ftanh 


?r 


<3tj 


) 


(31) 


Equation (24) is then used to transform from 4> to X, and T w (X) - can then be 
plotted along the impingement plate. The expression has been evaluated by numerical 
integration, and results are given in figure 3(a) for Pe = 5, 10, 20, and 50. The local 
Nusselt number is obtained from equation (30) and is plotted in figure 4(a). 


Simplified Solution Neglecting Axial Conduction 

In channel- flow heat-transfer analyses (ref. 11), the neglect of axial conduction 
produces only a small error for Pe between 10 and 100, and for Pe > 100 the error is 
negligible. Hence, it is worthwhile to examine here an approximate solution for which 
the axial conduction (conduction along the streamlines when transformed back into the 
jet geometry) is neglected. Compared with equation (31), this provides a much more 
simple solution that is convenient to evaluate and is accurate when Pe is sufficiently 
large. With this assumption the energy equation (25) simplifies to 
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(32) 


2 


0 * 


Pe 0T 
2 0 $ 


0 


which is of the same form as the transient-heat-conduction equation. The boundary 
condition (29) now applies along the axis $ = 0. The problem is the same as that for a 
plate of unit thickness that is initially at uniform temperature and that has one boundary 
kept insulated and one boundary with a heat input that varies with time. The solution can 
be obtained by using a superposition in time of the uniform -heat- flux solution in equa- 
tion (4) on page 112 of reference 6. The derivation is given in appendix B with the re- 
sult that the wall temperature is given by 


V*> ■ T 



Equation (24) is used to obtain 4> values corresponding to the desired X values, and 
the T W (X) are then evaluated by numerically integrating equation (33) for these 4>. The 
wall temperature distribution and corresponding Nusselt number distribution are plotted 
in figures 3(a) and 4(a) for various Pe. 


Simplified Solution Neglecting Axial Conduction and Upper Boundary 

The heat transfer in the vicinity of the stagnation region is sufficiently distant from 
the jet free streamlines that it is worth considering the solution of the energy equation 
where the effect of the upper boundary in figure 2 is neglected. As explained in the pre- 
vious section, the solution with axial conduction neglected is then the same as the solu- 
tion for the transient- heat- conduction equation - in this instance for a semi-infinite re- 
gion. From reference 6 (p. 76, eq. (9)) the solution for T w can be written as 

T (4>) - T = / f 

w ~ UW / U(*-if) 1/2 

Jo v 
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Inserting the U(<t>) for the jet configuration as given by equation (27) gives 


V*) - T 



( 34 ) 


The solution is thus a simplified form of equation (33) in which the infinite series ac- 
counted for the effect of the opposite wall. 

It is useful to evaluate equation (34) at the stagnation point as this is the region of 
highest heat transfer. As <3> — 0 in this equation, the 77 and $ - 7 7 will always be 
very small. For a small argument the hyperbolic tangent equals the argument, so equa- 
tion (34) becomes for small 



Let 77 = 77 /$ to obtain 


T '„<* - 0) - T„ 

W 77 



and then let 



to yield 


T w (*-0)-T„ =2 

77 





( 35 ) 
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From equation (30) 



Limiting Values for Large Values of Dimensionless Coordinate X 

For large X it is possible to obtain a simple limiting solution. From figure 1 it is 
evident that for large X the situation is a channel flow, and with viscous effects neglect- 
ed, the velocity is uniform. The channel is bounded at the top by the free streamline, 
which is assumed insulated, and at the bottom by the plate, which has a uniform heat 
addition. This region is treated as a channel flow in appendix C, and the result for the 
wall temperature at large X is 


V X) " T oo x+ ^ ( 37 ) 

This is compared with the solutions with and without axial conduction in figures 3(a) and 
4(a). 


DISCUSSION 

The condition considered here is for the impingement plate being heated uniformly, 
and the analysis yields the temperature distribution along the plate. The results are 
shown in figure 3; part (a) gives the temperature distribution along the plate, and 
part (b) gives the value at the stagnation point. The solid curves are the exact solution 
as given by equation (31), and as is typical of liquid metals, the results depend only on 
the single parameter Pe which does not involve viscosity. Numerical values are given 
in table I. As expected for a stagnation type flow, the highest heat transfer is at the 
stagnation point, and this is shown by the wall temperature being lowest at X = 0. The 
convective heat transfer should increase with the value of Pe, and the manner of this 
increase is revealed in figure 3. The heat- transfer characteristics can also be ex- 
pressed in the form of a local Nusselt number which is a reciprocal relation to the wall 
temperature, as shown by equation (30). The Nusselt number distribution along the 
plate and the stagnation Nusselt number are shown in figures 4(a) and (b). 
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An increase in convection, which corresponds to an increase in Peclet number, 
should diminish the effect of conduction along the streamlines, and the solution in equa- 
tion (33) should then apply. This relation is plotted as the long-dash curves in fig- 
ures 3(a) and 4(a). It is evident that the conduction term can be neglected, and the error 
will be within a few percent for Pe greater than about 10. Hence, unless Pe is some- 
what lower than 10, it is unnecessary to evaluate the exact solution as the simpler equa- 
tion (33) will suffice. 

At stagnation the neglect of the upper boundary and conduction along the streamlines 
leads to the following simple relations from equations (35) and (36): 

T w (X = 0)-T oo=2 J— (38a) 

K Pe 

Nu(X = 0) = <^5 (38b) 

As shown by the curves in figures 3(b) and 4(b), these relations are very good approxi- 
mations for Pe > 10. For low Pe the conduction along the streamlines at stagnation de- 
creases the heat transfer somewhat, as evidenced by the exact solution in figure 4(b) 
being below the solution neglecting heat conduction. Heat is conducted upstream into the 
flow as it approaches the plate, and this tends to thicken the thermal boundary layer 
near stagnation and decrease the heat transfer. Figure 3 shows that the wall tempera- 
ture is fairly uniform in the vicinity of stagnation. Hence, the stagnation results in 
equations (38) should also apply for the boundary condition of u ni form wall temperature. 

Equation (37) provides a simple asymptotic result for large X. The results show 
that this is a good approximation for X larger than about 6 with Pe greater than about 
10. At smaller X the asymptotic solution begins to deviate gradually from the exact 
results, and the approximation is not very good for X less than about 4 for the Pe 
values considered in this report. 

The analysis in reference 7 is concerned with gases and liquids, 0. 7 < Pr < 10, 
rather than the low Prandtl number range typical of liquid metals in the present report. 
The boundary condition at the impingement plate in reference 7 is a specified uniform 
wall temperature rather than a uniform heat flux like that treated in this report. For 
0 ^ Pr ^ 10 the viscous boundary layer is almost as thick or is thicker than the thermal 
layer, so that a boundary- layer analysis for the flow is required. The solution is ob- 
tained numerically with the jet potential flow solution used for the external flow. The 
analysis also includes the jet nozzle being at various distances from the plate. For 
large nozzle distances from the plate, as in the present report, the variations of Nusselt 
number along the plate have the same general trends as those obtained in this report. 
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1/2 0 373 

At the stagnation point the correlation obtained in reference 7 is Nu = 0. 51 Re ' Pr 

as compared with Nu = 0. 707 Pe^ = 0. 707(RePr)* // ^ for the present report. The cor- 

1 /2 0 1 27 

relation in reference 7 can be rearranged into Nu = 0. 51(RePr) ' /Pr = '* 

1/2 0 127 

0. 51 Pe ' /Pr . Viscosity appears only in the Prandtl number, and since the 
Prandtl number is to a small power, the heat-transfer behavior at the stagnation point 
is not very viscosity dependent. 


CONCLUSIONS 

An analysis has been performed combining free jet theory and a transformation of 
the energy equation into the potential plane to obtain the heat-transfer behavior of an 
impinging liquid-metal jet. The jet is two-dimensional in Cartesian coordinates and 
strikes a plate that is uniformly heated. The temperature distribution and heat-transfer 
coefficient are evaluated along the plate as a function of the Peclet number Pe. As 
would be expected, the maximum heat-transfer coefficient is at the stagnation point. At 
a distance larger than about 6 jet half-widths from stagnation and for Pe > 10, the re- 
sults can be approximated quite well by an asymptotic relation that considers the jet 
region to act like a channel flow. Heat conduction along the streamlines is not very sig- 
nificant when Pe is greater than 10. For this Pe range, the Nusselt number Nu at stag- 
nation can be approximated quite well by Nu = yj Pe/2 . This result should also apply for 
the boundary condition of uniform wall temperature because the present solution yielded 
almost uniform wall temperatures in the vicinity of stagnation. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 7, 1973, 

501-24. 
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APPENDIX A 


SOLUTION OF ENERGY EQUATION 

The solution can be obtained by starting with the equation for a stationary point 
source of heat on one side of a moving slab with the other side of the slab insulated. By 
superposition the point source is used to obtain the solution for a line source. Then a 
superposition of line sources can be used to build up the heat source distribution speci- 
fied by the boundary condition equation (27) along the positive axis. 

For a slab of unit thickness, 0 < 4* <1, moving with a unit velocity (fig. 2) the tem- 
perature distribution produced by a unit point heat source fixed at the origin, 4> = 4' = 
|=0 (the £ is a coordinate axis perpendicular to $ and 'k), is obtained from refer- 
ence 6, page 268, equation (7), as 



Integrating over £ yields the temperature distribution for a line source along the £ 
axis 



To carry out this integration the result from reference 12, page 417, is used that 
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(A3) 



2 M r(ju + i) 

aM+y-M- 1 




a > 9 


o 

Since the integrand in equation (A3) contains only ^ , the result should be the same for 
negative j as for positive Hence, to allow for negative f the result in equa- 
tion (A3) should more properly be written as 




2 M r (] a + 1 ) 




(A4) 


A specific case for v - 0 and /i = -1/2 is 



Using the relations = K i/2 J r (V2) = Vir, and V(ff/2^) Kqy 2 (?) = ( n /2f) e - ? 

(see ref. 13, p. 444) yields the simplified form 



(A5) 


From symmetry equation (A2) can be written with integration limits 0 — °° as 
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Equation (A6) is integrated by use of equation (A5) to yield 


r 





TT 


1 J_ e - (Pe/4) l*l 

2 Pe 




CO 



m=l 


■ (Pe/4) [l+(4m zr/Pe) 2 ] 1/2 1 $ 


n e 


2 





cos(mn'\I r )y e 


$Pe/4 


Simplifying and evaluating at the wall ^ = 0 give 


V*> • T 


/ 

_ _2_ e (Pe/4)(*-|*|) 
Pe 




oo 



m=l 


2 

(Pe/4) j$- | $| |l+(4m?r/Pe) 2 ] | 

e ^ 


1/2 

1 /4m7r\ 

y 

\ Pe / 


(A7) 


Equation (A7) gives the wall-temperature distribution produced by a line source of 
unit strength per unit length along the £ axis in figure 2. Superposition can now be 
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applied along the positive $ axis to obtain the effect of the source distribution 1/U(<£) 
which resulted from the imposed uniform heat flux in the physical plane being trans- 
formed into the potential plane. If H(<3>) is the response to a unit line source, and F($) 
is the distribution of sources, the wall- temperature distribution can be obtained by 
superposition as 







dr? 


(A9) 


The U(tj) is inserted from equation (27) and the integral is divided into two regions de- 
pending on whether r\ < <£ or rj > $ to obtain the final form 
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- (Pe/4)(7j - <£) jl+ |l+ (4mff/Pe) 2 J 

a v 4 


^tanh!5 



\ 

I77 

I 


(A10) 
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APPENDIX B 


SOLUTION NEGLECTING AXIAL CONDUCTION 


In reference 6, page 112, the transient temperature solution is given for a heat flux 
suddenly imposed at one surface of a slab with the other surface of the slab insulated. 
Since equation (32) has the same form as the transient heat conduction equation, this 
solution can be applied after changing the nomenclature. The <f> is analogous to time in 
the transient problem, and since the heat flux varies with $ in the present case, the 
solution is the same as a transient solution in which the imposed heat flux varies with 
time. The superposition theorem for the variable- heat- flux solution is 


T w (*) - T 


00 



1 dF 
U($ - T]) dr] 


dr] 


(Bl) 


where l/U($) is given by equation (27), and F is the response to a suddenly imposed 
heat flux of unit magnitude. 

From the solution in reference 6 the F is given by 


OO 



m=0 


r 




< ierf c 

s. 


m 



+ ierfc 


m + 1 



where 


ierfc 4 = — e 


-i 1 


- 4 erfc 4 


The derivative of this function is given by 

A ierfc | = dgerfcjj dg_ = ( _ erfc ^ d£ 
d$ d4 d$ d$ 
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Using this relation gives 


oo 



When the definitions of ierfc 4 and erfc 4 are used, equation (B2) reduces to 


dF = ( 2 \ 1/2 J_ V [ -m 2 /(2/Pe)4> -(m+l) 2 /(2/Pe)$] /- 

d* UPe/ -I 

m=0 

Substitute equation (B3) into equation (Bl) and use equation (27) for U to obtain the de- 
sired final result for the wall-temperature distribution 



(B4) 
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APPENDIX C 

LIMITING HEAT TRANSFER AT LARGE VALUES 
OF DIMENSIONLESS COORDINATE X 

For large x the flow as shown in figure 5 is a flow with uniform velocity in a chan- 
nel with one wall heated and one wall insulated. An overall heat balance is now taken on 
the flow region from the fluid entrance at large y to the region at large x, where the 
velocity is uniform and equal to v^ . This yields for large x 

V s =bv 00 pC p [t(x) - tj (Cl) 


where t (x) is the mean temperature for uniform velocity (see eq. (C7)) obtained as 

- 1 r b - 

t (x) = ^ Jq t(x, y)dy. Equation (Cl) determines the t (x) since q w is a specified quan- 
tity. 

A heat balance on a differential element in the fluid gives 


k 




' v “ pc pf x 


(C2) 


For fully developed heat transfer with uniform heat addition the temperature distribution 
rises linearly with x and hence has the form 


t(x, y) = CjX + f(y) (C3) 

2 2 

The axial conduction term k(3 t/3x ) is thus zero, and substitution of equation (C3) into 
(C2) yields 


j2r 

k— =v 00 pCC 1 (C4) 

dy 2 


Since in the fully developed region the shape of the temperature distribution is not chang- 
ing, equation (Cl) yields 


3t _ 3t ^w 
3x 3x bv^pCp 


(C5) 
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Then combining equations (C5) and (C4) gives 


1 ? 

dy 


d 2 f = % 
2 kb 


Integrate once and apply the boundary condition that 3t/ dy = 3f/3y = 0 at y = b to ob- 
tain 


df 


=— (y-b) 
dy kb 


A second integration results in 


f = — (^-by] + -C 2 
kb \ 2 / * 


The integrated average temperature is 

fo ut d y v co Jo 


t (x) = : 


b 

o td y 


f b f b b 

Jo ud y v =o Jo d y 


_i f b 

~7 JO 


t dy 


fO “ v oo Jo 

Substituting for t by use of equations (C3), (C5), and (C6) gives 


t(x) 

b 


i{ 


V 


bv^pC 


t(x) - - 1 


q / 

' 2 


y 

kb \ 

,2 

_ b^ 

kb 

3 


dy 


(C6) 


(C7) 


Using t(x) from equation (Cl) gives 

V 


bv pC 

P 


ax a b 
b Voo pC 3k 
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so that 


C 2 " t oo + 


q b 
3k 


(C8) 


Combining equations (C3), (C5), (C6), and (C8) gives the temperature distribution 



V 

kb ' 

U 


bV coP C p 

Evaluating at y = 0 to obtain the wall temperature yields 


q b 

4 w 

3k 


V x ) " ^ = 

In dimensionless form the final result is 

T w (X) - T 


V , q w b 

bV ooP C p 3k 



(C9) 
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Dimensionless wall temperature, T W (XI - 




Dimensionless distance from stagnation, X 



(a) Distribution along impingement plate. 


(b) Dimensionless wall temperature at stagnation (X = 0); 
can be approximated with very small error for Pe > 10 


by T w -Tcc= 2V2JP5. 


Figure 3. - Dimensionless wall temperature for impinging jet. 
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Local Nusselt number, Nu(X) 


Solution Equation 


Exact (31) 

Without conduction (33) 

along streamlines 

Asymptotic for large X (37) 




(a) Local Nusselt number along impingement plate. 


(b) Nusselt number at stagnation (X * 0); can be approxi- 
mated with very small error for Pe>10 by Nu ■ 
(Pe/2) 1/2 . 


Figure 4. - Nusselt number for impinging jet. 



Figure 5. - Geometry for derivation of heat-transfer behav- 
ior at large x. 
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